Energy and enstrophy fluxes in the double cascade of 2d turbulence 
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High resolution direct numerical simulations of two-dimensional turbulence in stationary condi- 
tions are presented. The development of an energy-enstrophy double cascade is studied and found 
to be compatible with the classical Kraichnan theory in the limit of extended inertial ranges. The 
analysis of the joint distribution of energy and enstrophy fluxes in physical space reveals a small 
value of cross correlation. This result supports many experimental and numerical studies where 
only one cascade is generated. 



PACS numbers: 

The existence of two quadratic inviscid invariants is 
the most distinguishing feature of Navier Stokes equa- 
tions in two dimensions. On this basis, in a remarkable 
paper in 1967 [1], Kraichnan predicted the double cas- 
cade scenario for two-dimensional turbulence: an inverse 
cascade of kinetic energy E — l/2(u 2 ) to large scales and 
a direct cascade of enstrophy Z = l/2(w 2 } to small scales 
(u> = V x v is the vorticity). In statistically stationary 
conditions, when the turbulent flow is sustained by an 
external forcing acting on a typical scale if a double cas- 
cade develops. According to the Kraichnan theory, at 
large scales i.e. wavenumbers k <C kf ~ ij 1 , the energy 
spectrum has the form E(k) ~ e 2 / 3 fc -5 / 3 , while at small 
scales, k 3> kf, the prediction is E(k) ~ r] 2 / 3 k~ 3 , with a 
possible logarithmic correction [2]. Here e and r\ ~ fc 2 £ 
are respectively the energy and the enstrophy injection 
rate. 

Despite the importance of two-dimensional turbulence 
as a model for many physical flows [3, 4] and, more in gen- 
eral, for non-equilibrium statistical systems [5], a clear 
evidence of the coexisting two cascades on an extended 
range of scales is still lacking. The inverse energy cas- 
cade has observed in many laboratory experiments [6] 
and in numerical simulations [7, 8, 9] with a statistical 
accuracy which has revealed the absence of intermittency 
corrections to the dimensional scaling [10]. For what con- 
cerns the direct cascade, earlier numerical simulations 
and experiments report spectra slightly steeper than fc -3 
[11, 12], while more recent investigations at high resolu- 
tion are closer to Kraichnan prediction [13, 14, 15, 16]. 
It is important to remark that in presence of a large scale 
drag force (always present in experiments and sometimes 
also in numerics) one indeed expects a correction to the 
classical exponent —3 [17, 18]. 

Two recent experimental papers have been devoted to 
the study of the double cascade [19, 20]. Their results 
are substantially consistent with the classical scenario of 
Kraichnan, although the extension of the inertial range 
(in particular for the inverse cascade) is limited. Here we 
present high resolution (up to 16384 2 ) direct numerical 
simulations of forced 2D Navier-Stokes equations which 



reproduce with good accuracy both the cascades simulta- 
neously. Most of the injected energy flows to large scales 
(where it is removed by friction dumping) while enstro- 
phy cascades to small scales (there removed by viscosity). 
We find strong numerical indications that the classical 
Kraichnan scenario is recovered in the limit of two ex- 
tended inertial ranges, although we are unable to rule 
our the possibility of small corrections in the direct cas- 
cade. By looking at the two fluxes in physical space, 
we find a relatively small value of the cross correlation 
among them. This result is interpreted in favor of the 
possibility of generating a single cascade, independently 
on the presence of the second inertial range. 

The 2D Navier-Stokes equation for the vorticity field 

is 

d t ui + v ■ Vw = v Auj — auj + Af, (1) 

where v is the kinematic viscosity and a is a linear friction 
coefficient (representing bottom friction or air friction) 
necessary to obtain a stationary state. The forcing term 
/ is assumed to be short correlated in time (in order to 
control the injection rates) and narrow banded in space. 
Specifically, we use a Gaussian forcing with correlation 
function (/(r,i)/(0,0)) = 5{t)Fl)exp{-{r/i s f/2) in 
most of the simulations. For the simulations at resolu- 
tion 16384 we use a different forcing restricted to a nar- 
row shell of wavenumbers. Numerical integration of (1) 
is performed by pseudo-spectral, fully dealiased, parallel 
code on a doubly periodic square domain at resolution 
iV 2 . Statistical quantities are computing by averaging 
over several large-scale eddy turnover times in stationary 
conditions (only over a fraction of eddy turnover time for 
the 16384 run because of limited resources). Table I re- 
ports the most important parameters of our simulations. 

One of the simplest information which can be obtained 
from Table I is related to the energy-enstrophy balance. 
At N — 2048 only about one half of the energy injected is 
transfered to large scales where it is removed by friction 
at a rate e a — 2aE. This fraction increases with the res- 
olution and becomes about 95 % for the N — 16386 run. 
The remaining energy injected is dissipated by viscosity 
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N v_ a Ijjld Ea/ei £ v /ei rja/rji r/ v /rji 

2048 2 x 1(T 5 0.015 26.2 0.54 0.46 0.03 0.97 

4096 5 x 10~ 6 0.024 52.3 0.82 0.18 0.08 0.92 

8192 2 x 10~ 6 0.025 80.5 0.92 0.08 0.10 0.90 

16384 1 x 10~ 6 0.03 114.2 0.95 0.05 0.12 0.88 

TABLE I: Parameters of the simulations. N spatial resolu- 
tion, v viscosity, a friction, L box size, t$ = i/100 forcing 
scale, Id — v 1,/2 /ri^ 6 enstrophy dissipative scale, si energy 
injection rate, e„ viscous energy dissipation rate, e a friction 
energy dissipation rate, rji enstrophy injection rate, i) v viscous 
enstrophy dissipation rate, rj a friction enstrophy dissipation 
rate. 

at scales comparable with the forcing scale and at a rate 
proportional to v (which thus decreases by increasing the 
resolution) . 

On the other side, most of the enstrophy (around 90 %) 
follows the direct cascade to small scales, where it is dis- 
sipated by viscosity. We observe a moderate increase of 
the large scale enstrophy dissipation r\ a increasing the 
resolution: this is a finite-size effect due the increase of a 
with TV (see Table I) necessary to keep the friction scale 
£ a ~ el/ 2 a~ 3 / 2 constant with increasing e a . 

In Figure 1 we plot the fluxes of energy and enstro- 
phy in wavenumber space. Observe that because change 
the resolution keeping the ratio L/lf constant, the only 
effect of reducing the resolution on the inverse cascade 
is the decrease of the energy transfered to large scales 
(being £„=£/ — e v with e v proportional to v) while the 
extension of the inertial range is almost unaffected. The 
qualitative difference at k ~ kf for the run at N = 16384 
is due to the different forcing implemented in this case, 
while statistical fluctuations are a consequence of the 
short time statistics. These results confirm the robust- 
ness of the energy inertial range regardless of the viscous 
dissipative scale, a further justification of many simu- 
lations of the inverse cascade in which, because of the 
limited resolution, the forcing scale is very close to the 
dissipative scale. 

At variance with the inverse cascade, the direct en- 
strophy cascade is strongly affected by finite resolution 
effects. This is not a surprise because, by keeping if 
fixed, the extension of the direct cascade is simply pro- 
portional to AT. As it is shown in Figure 1, we observe a 
range of wavenumbers with almost constant flux Hz(k) 
only for the runs at N > 8192. 

Figure 2 shows the energy spectra computed for the dif- 
ferent runs. We remark again that the only effect of finite 
resolution on the inverse cascade is the reduction of the 
energy transfered to large scales, while the Kolmogorov 
scaling fc~ 5 / 3 is always observed with a Kolmogorov con- 
stant C ~ 6 [10] virtually independent on resolution. 
The effect of finite resolution is, of course, more dra- 
matic on the enstrophy cascade range. We observe here 
a significative correction to the Kraichnan spectrum fc~ 3 
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Fourier space at resolutions 2048, 4096, 8192 and 16384. Res- 
olution increases from smaller to larger absolute values of the 
minimum/maximum. At resolution 16384 statistics is com- 
puted on a single frame. 

even for the 16384 run, where measure a scaling expo- 
nents close to —3.6. We remark a similar steepening 
of the spectrum has been observed even for simulations 
with a more resolved direct cascade range (here we have 
kmax/kf ~ 55 at the highest resolution). Despite these 
difficulties, there is a clear indication that the correction 
to the exponent is a finite-size effect which eventually dis- 
appears by increasing the extension of the inertial range 
(see inset of Fig. 2). The conclusion of these consider- 
ations is that a A; -3 spectrum in stationary solutions of 
(1) could be achieved only by taking simultaneously the 
limits L/£f — > oo and if /id — ► °o (i.e. for vanishing a 
and v). 

A better understanding of the physical mechanism of 
the cascades can be obtained by looking at the distri- 
bution of fluxes in space. This can be obtained by us- 
ing a filtering procedure recently introduced and applied 
separately to the direct [21] and to the inverse [22] cas- 
cades. Thanks to the resolution of the present simula- 
tions, we are able to analyze jointly both the cascade 
and the correlations among them. Following [21], we 
introduce a large-scale vorticity field uj r = G r * u> ob- 
tained from the convolution of u> with a Gaussian filter 
G r , and a large-scale velocity field v r = G r * v. From 
these definitions, balance equations for the large-scale en- 
ergy e r (x,t) = l/2|t; r | 2 and enstrophy z r (x,t) — 1/2lu% 
densities are easily written (with a compact notation): 

d t (e r , z r ) + V • J^ r) = -Ii^ r) - + (2) 

where Jr 6 ' 2 -* represent the transport of large-scale en- 
ergy and enstrophy, D'f'^ and Fr e ' r ^ represent the 
large-scale dumping and forcing. The energy /enstrophy 
fluxes ni 6 ' 2 ' (x, t), representing the local transfer of en- 
ergy /enstrophy from large scale to scales smaller than r, 
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FIG. 2: Energy spectra for the two simulations for the dif- 
ferent resolutions (resolution increases with the extension of 
the tail). Dashed and dotted lines represent the two predic- 
tions Ck~ r '/ 3 with C = 6 and fc~ 3 respectively. Inset: cor- 
rection to the Kraichnan exponent —3 measured in the range 
100 < k < 400 as a function of viscosity. 

are given by: 

rif\x,t) = -(r aP ) r V a (vp)r (3) 
II^OM) = -{a a ) r V a uo r (4) 

where (r Q/3 ) r = {v a vp) r - {v a ) r (vp) r and (a a ) r = 

(v a Uj) r - {v a ) r U) r . 
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FIG. 3: Average energy (lower, continuous) and enstrophy 
(upper, dotted) fluxes in physical space at resolution N = 
8192. The energy flux is multiplied by a factor k 2 f = 10 4 for 
better visualization. 

Fluxes (3-4) are expected to have a non-zero spatial 
mean in the inertial range of scales of irreversible turbu- 
lent cascades: in particular a mean negative energy flux 
for r > if (inverse cascade) and a mean positive enstro- 
phy flux for r < if (direct cascade). Figure 3 shows the 
averaged physical fluxes (ni e ' z ^) as functions of the scale 



r. The two cascades are evident, although the range of 
constant flux is reduced with respect to the spectral case 
in both the cascades (see Fig. 1). 

Local fluxes are strongly inhomogencous in physical 
space: there are relatively small regions of intense (posi- 
tive and negative) flux in both the energy and enstrophy 
inertial ranges. Figure 4 shows two snapshots of the en- 
ergy and enstrophy fluxes, computed from the same field 
at two different scales n and r-i corresponding to the min- 
imum of energy flux and the maximum of enstrophy flux 
respectively (see Fig. 3). An interesting information, ob- 
tained from Fig. 4 at a qualitative level, is that the most 
intense energy and enstrophy fluxes appear on different 
physical region without any apparent correlation. 




FIG. 4: Snapshot of energy (left) and enstrophy (right) fluxes 
from the same vorticity field. Energy flux (3) is computed 
at scale ri = 0.0251/ and enstrophy flux at ri = 0.0025L, 
roughly corresponding to the minimum and maximum in 
Fig. 3. White and black correspond to positive and negative 
values respectively. 

Figure 4 shows that both positive and negative fluxes 
are observed: locally both energy and enstrophy can go 
to smaller or larger scales. Figure 5 shows the probability 
density functions of the two fluxes. As observed in previ- 
ous studies [21, 22], the shapes of pdf's are nearly sym- 
metric, confirming the qualitative picture inferred from 
Fig. 4. The mean values in both cases is the result of 
strong cancellations as its ratio with the standard devia- 
tion is —0.22 for the energy flux and 0.16 for the enstro- 
phy flux. The skewness is also small as it is about 0.9 
and 1.2 for the energy and enstrophy fluxes respectively. 

Figure 6 shows the joint probability density function 
ptn^.n^) for the same scales r\ and r^- This pdf 
is not far from the product of the marginal distribu- 
tions shown in Fig. 5, a condition for independence. In- 
deed, the cross correlation among IT^ and Ilif is only 
C( r ii r 2) — —0.15. Of course, it is very different if we 
consider the correlation at the same scale, for which we 
find C(r,r) ~ 1 despite the fact that one of the fluxes 
has zero mean. A possible interpretation of the observed 
small value of correlation is in favor of the classical pic- 
ture of independence of the two cascades which is here 
obtained at a local level. Therefore, our result is, a pos- 
teriori, a support to many experimental and numerical 
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FIG. 5: Probability density function of energy nil' ( left ) and 
enstrophy 11^ (right) fluxes normalized with their standard 
deviations. Resolution is N = 8192, filtering scales are r\ = 
0.025L and r 2 = 0.0025L. 

studies in which, due to finite size constraints, only one 
cascade is realized. 
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FIG. 6: Joint probability density function pfTli^ , n^') of en- 
ergy (x-axis) and enstrophy (y-axis) fluxes. Scales n and r-2 
as in Fig. 5. Contours are plotted on a logarithmic scale while 
fluxes are normalized with their standard deviations. 

In conclusion, we have presented statistical analy- 
sis of high resolution direct numerical simulations of 
2D Navier-Stokes equations which clearly reproduce, for 
the first time, the double cascade scenario predicted by 
Kraichnan almost 40 years ago. 

Simulations were performed on the IBM-CLX cluster 
of Cineca (Bologna, Italy) and on the Turbofarm cluster 
at the INFN computing center in Torino. 
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